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The quantum corrections to the energies of the T point optical phonon modes (Kohn anomalies) in 
graphene nanoribbons are investigated. We show theoretically that the longitudinal optical modes 
undergo a Kohn anomaly effect, while the transverse optical modes do not. In relation to Raman 
spectroscopy, we show that the longitudinal modes are not Raman active near the zigzag edge, while 
the transverse optical modes are not Raman active near the armchair edge. These results are useful 
for identifying the orientation of the edge of graphene nanoribbons by G band Raman spectroscopy, 
as is demonstrated experimentally. The differences in the Kohn anomalies for nanoribbons and for 
metallic single wall nanotubes are pointed out, and our results are explained in terms of pseudospin 
effects. 



I. INTRODUCTION 

Graphene nanoribbons (NRs) are rectangular sheets 
of graphene with lengths up to several micrometers 
and widths as small as nanometers NRs can be 
regarded as unrolled single wall nanotubes (SWNTs). 
Since SWNTs exhibit cither metallic or semiconducting 
behavior depending on the diameter and chirality of the 
hexagonal carbon lattice of the tube^ it is expected that 
the electronic properties of NRs depend on the width and 
"chirality" of the edge.— In fact it has been predicted 
that their electronic properties near the zigzag edge are 
quite different from those near the armchair edge . 8 ' 9 ' 10 ! 11 
Thus the characterization of the NRs as well as SWNTs 
is a matter of prime importance. 

Raman spectroscopy has been widely used for the 
characterization of SWNTs 1 ^^^ 1 ^! 8 .^^* 21 . and 
graphenes . 22 ' 23 ' 24 ' 25 ' 26 ' 27 Recently, it has been shown 
that the frequencies and spectral widths of the T point 
optical phonons (called the G band in Raman spectra) 
depend on the position of the Fermi energy E-p and the 
chirality of the metallic SWNT.^ 2 ^ 3 ^ The Fermi en- 
ergy dependence of the Raman spectra can be used to 



determine the position of the Fermi energy, and the chi- 
rality dependence of the Raman spectra provides detailed 
information on the electronic properties near the Fermi 
energy of metallic SWNTs. These dependences originate 
from the fact that the conduction electrons of a metal 
partly screen the electronic field of the ionic lattice. Kohn 
pointed out that the ability of the electrons to screen the 
ionic electric field depends strongly on the geometry of 
the Fermi surface, and this screening leads to a change in 
the frequency of a specific phonon and an increase in its 
dissipation (the Kohn anomalies) . 32 The Kohn anomalies 
(KAs) in graphene systems are unique in the sense that 
they can occur with respect to the T point phonons while 
the KA in a normal metal occurs with respect to phonons 
with 2kp where hp is the Fermi wave vector. The unique- 
ness of graphene comes from the geometry of the Fermi 
surface given by the Dirac cone ^i34^^ML^i^iM. 

Since the geometry of the Fermi surface of NRs, and 
the energy band structure of NRs depend on the orien- 
tation of the edge^ ' 10 ' 11 one may expect that the KAs 
of NRs depend on the "chirality" of the edge like the 
KAs of metallic SWNTs. In this paper, we study KAs 
in graphene NRs with zigzag and armchair edges. A NR 



with a zigzag (armchair) edge is hereafter referred to as 
a Z-NR (an A-NR) for simplicity (see Fig. QJa) for an N 
Z-NR with a width W = N£ where £ = 3a cc /2 and a cc 
(= 0.142nm) is the bond length). We will show that the 
transverse optical (TO) phonon modes do not undergo 
KAs in both A-NRs and Z-NRs. The dissipation of the 
longitudinal optical (LO) phonon modes in Z-NRs is sup- 
pressed as compared to those in A-NRs. We also show 
that the LO (TO) modes are not Raman active in Z-NRs 
(A-NRs), and that the KAs should be observed in only 
A-NRs. This fact is useful in identifying the orientation 
of the edge of NRs by G band Raman spectroscopy. 

It is noted that the D band which consists of an in- 
tervalley K point phonon has been used to characterize 
the orientation of the edge of graphite and graphene. 
Pimenta et al^- observed that the intensity of the D 
band near the armchair edge of highly ordered pyrolytic 
graphite (HOPG) is much stronger than that near the 
zigzag edge of HOPG. This is confirmed for single layer 
graphene by the experiments of You et al— who also 
show that the D band has a very strong laser polarization 
dependence. However, a strong D band intensity appears 
only ~ 20 nm from the edges, so that the observation of 
this effect requires a precise experimental technique^ 

This paper is organized as follows. In Sec. Hit we study 
KAs in Z-NRs and A-NRs. In Sec. DTD we point out that 
the LO modes in Z-NRs (TO modes in A-NRs) are not 
Raman active and show experimental results. In Sec. \TV\ 
the mechanism of the edge dependent KAs and Raman 
intensities is explained in terms of the pseudospin. A 
discussion and summary are given in Sec. IVland Sec. IVII 



II. KOHN ANOMALIES 

A. Zigzag NRs 

KAs are relevant to the electron-phonon (el-ph) matrix 
element for electron-hole pair creation. The electron-hole 
pair creation should be a vertical transition for the T 
point optical phonon, and the KA effect, such as an in- 
crease in the dissipation of the phonon, is possible only 
when the energy band gap of a NR is smaller than the 
energy of the phonon (about 0.2eV— ). First, we study 
the energy band structures of Z-NRs whose lattice and 
phonon modes are shown in Figs.[lja) and (c,d), respec- 
tively. Z-NRs have a metallic energy band structure re- 
gardless of their widths as shown in Fig. [Ha) . The metal- 
licity of Z-NRs is due to the edge states forming a fiat 
energy band at E = 0. 9 Similarly, armchair SWNTs have 
a metallic energy band regardless of their diameters as 
shown in Fig. [^b)^ It is interesting to imagine that an 
N = 2n— 1 Z-NR can be obtained from an (n, n) armchair 
SWNT by cutting the circumferential C-C bonds along 
the tube axis as shown in Fig. [T](b) . For example, from 
a (5, 5) armchair SWNT, we get an N = 9 Z-NR. The 
metallicity of armchair SWNTs is preserved in Z-NRs by 
the edge states4£ 
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FIG. 1: (a) The lattice structure of a Z-NR. The solid 
(empty) circles denote the A (B) sublattice. We use integers 
I G [0, M] and J G [0, N] for the axes. The width (length) W 
(L) of a Z-NR is given by N£ (oM where a = V3a cc ). (b) An 
(n, n) armchair SWNT is cut along its axis and flattened out 
to make an N = 2n — 1 Z-NR. There is a metallic energy band 
in both structures (see Fig. [2ja) and (b)). (c,d) The phonon 
eigenvectors of the F point LO and TO modes are illustrated. 
The LO mode satisfies 570,1 = and 570,2 = — 570,3. The TO 
mode is characterized by 570,1 = —2570,2 and 570,2 = 570,3, 
and 570,1, 570,2, 570,3 are defined by the inset. 



® JV=9Z-NR CO (5,5) SWNT ( c ) 




FIG. 2: The energy band structure of an N = 9 Z-NR (a) and 
that of a (5, 5) armchair SWNT (b). (c) |Vl/<57| as a function 
of ka. The matrix elements of the vertical transitions, denoted 
by the arrows (1,2) in (a) and (b), are shown as the solid and 
dashed curves, respectively. 



The phonon eigenvector of the LO (TO) mode is par- 
allel (perpendicular) to the zigzag edge as shown in 
Fig. [TJ;(d). By the displacement of a C-atom, a bond 
length increases or decreases depending on the position 
of the bond. A change of bond length causes a change 
of the three nearest-neighbor hopping integrals from —70 
(= 2.7eV) to —70 + 57o,a (a = 1,2,3) (see the inset to 
Fig. []Tc,d)). The electron eigen function in the presence 
of 070, a is given by a linear combination of those in the 
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absence of <57o, a . In other words, a shift <57o, a works as 
a perturbation by which an electron in the valence band 
may be transferred to a state in the conduction band. 
This is an electron-hole pair creation process due to a 
lattice deformation. We derive the el-ph matrix element 
as follows. 

Let V'a' 7 (V'b J ) denote the wave function of an electron 
at an A-site (B-site) at site 13 {I 1 J where V = I + 1/2). 

Then X)/(V'a' / )*^7o,3V'b J 1S the amplitude for the pro- 
cess that an electron at the B-sites with I'J is trans- 
ferred by the perturbation Sjo, 3 into the A-sites with IJ 
(see Fig.[IJa)). Note that (V'aO* indicates complex con- 
jugation of ?/> A J and Sjo,3 for the LO and TO modes is 
constant. By introducing the Bloch function (</>a>%) 
as Vl J = {e lI(ka) /VM)4>i and V B ' J = {e lI ' [ka) /VM)^ 
where k is the wave vector along the zigzag edge, we 
obtain £j(^a )**70,3< J = «)^ 7 o,3e lfca/2 0B- % d e- 
riving the matrix elements for $70. 2 and <5 7 o,i in a similar 
manner, we obtain the el-ph matrix element for a vertical 
electron-hole pair creation process such as V + U, where 
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l\ ( e 1 570,2 + e l ^670,3^ (&A 
c. 
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(1) 



V and U represent the el-ph interaction for the £70. 2 
(^7o, 3) perturbation which acts for the same J and the 
^70,1 perturbation which acts for the nearest neighbor 
pair of J and J', respectively. For V in Eq. {!]), c.c. repre- 
sents the complex conjugation of e~ l ~ ^70,2 + e I_ 2~<57o,3. 
It is also noted that the minus signs in front of cAg in 
Eq. {1} come from the fact that the Bloch function with 
energy E in the conduction band is given by (</> A ,— 0g) 
when the Bloch function with energy — E in the valence 
band is (</> a ,<^>b). This is a property of the nearest- 
neighbor tight-binding Hamiltonian with two sublatticcs. 
The matrix elements V and U in Eq. |T]) can be rewritten 



<^7o,2 = — ^7o,3 for the LO mode, while c>7o,i = — 2^70,2 
and 570 2 = ^7o 3 for the TO mode (see Fig. QJc) and 
(d)). 

Thus, for the LO mode, both Vt and U vanish because 
the LO mode satisfies £70, 2 + ^70,3 = and c>7o,i = 0, 
respectively. The non-vanishing matrix element for the 
LO mode is given by Vl only. By setting $70, 2 = — £70,3 
and introducing a shift £7 due to a bond stretching as 
£70,3 = <57COs(7r/6) in Eq. ([3]), we have 



V L = 2V3i<5 7 sin ( y 



J 



(5) 



V = V t + Vl 



(2) 



From Eq. ([5]) , it is understood that the electron-hole pairs 
around ka = are not excited since \Vl\ is proportional 
to sin(fca/2). Moreover, since the wave function of the 
the edge states appears only on one of the two sublat- 
tice in the hexagonal unit cellj2 we have J2j 'Pa^b ~ 
in Eq. ([5]) for the edge states. This suppresses the el- 
ph matrix element of electron-hole pair creation for the 
edge states. These facts can be checked by a numerical 
calculation as shown in Fig. [2]Jc) where we plot |Vl/#7| 
as a function of ka by the solid curves for the two low- 
est energy electron-hole pair creation processes which are 
denoted by the arrows in Fig. [2a). In Fig. [2jc), we 
see that the edge states appearing as a flat energy band 
at 27r/3 < ka < it do not contribute to electron-hole 
pair creation (solid line 1). In Fig. [DJc), we also plot 
|Vl/^7| for a (5,5) armchair SWNT for comparison by 
the dashed curves for the two lowest energy electron-hole 
pair creation processes which are denoted by the arrows 
in Fig. HJb). As for the lowest energy electron- hole pairs 
(the dashed line 1), |Vl| increases with increasing ka due 
to sin(fca/2) in Eq. ([5|). This indicates a constant value of 
12 j 4>a4>b fo r the case of armchair SWNTs. In Fig. EJc), 
we see that the matrix element of the next lowest energy 
electron- hole pairs vanishes at kg satisfying de/dk\k a = 
(van Hove singularity) for both the Z-NR and armchair 
SWNT. The same behavior is observed for higher sub- 
bands, and a large density of states due to the van Hove 
singularities of the sub-bands is not effective in producing 
the electron- hole pair. 

For the TO mode, Vl vanishes owing to $70,2 — ^7o,3 = 
0. Moreover, it can be shown that 



where 



V T = 2i (<5 7 o,3 + £70,2) cos ( ) £ Im [0 A *«^] 



Vl = 2i (570,3 - 6^,2) sin (^j £ Re [<&4&] , 



and 



(3) 



(4) 



By assuming that the perturbation <57o, a is proportional 
to a change of the bond length, we have £70,1 = and 



Im 



f>i*<f>i] = 0, In# A +1 *0 B ] = 0, 



(0) 



since an analytic solution of 
in Ref. 13 as 



>b) for Z-NRs is given 



— sin Jcf> 
9 



sin(J + 1) 



CM), 



9 



sin(J+ 1) 



(7) 



Here g = 2 cos(fca/2), and C(g, </>) is a normalization con- 
stant, <f> is the wave number which is determined by the 
boundary condition: A +1 = 0, and e(g, cf>) is the en- 
ergy eigenvalue in units of —70- The energy dispersion 
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relation is given by e(g,cj)) 2 — g 2 + 1 + 2gcos</>. 44 Since 
C* {g , 4>)C (g , (j>) is a real number, we get Eq. ©. As a 
consequence, we have Vt = and U = in Eqs. © and 
(0D. Thus, the TO modes give both V = and U = for 
any vertical electron-hole pair creation process, and the 
TO modes decouple from the electron-hole pairs. This 
shows the absence of the KA for the TO modes in Z- 
NRs. 

A renormalized phonon energy is written as a sum 
of the unrenormalized energy Uujq and the self-energy 
n(a>o)- Throughout this paper, we assume a constant 
value for Tiloq as Tiloq = 1600cm -1 (0.2eV) both for the 
LO and TO modes. The self-energy is given by time- 
dependent second-order perturbation theory as 



IVU 5 



hu - E eil + iT/2 Tiui + E c h + iT/2 



X (A - /c) ■ 



(8) 



where the factor 2 comes from spin degeneracy, /h. c — 
(1 + exp(/3(i5h,c — Epj)^ 1 is the Fermi distribution func- 
tion, E e (£?h) is the energy of an electron (a hole), and 
i^ch = E c — Eh is the energy of an electron-hole pair. In 
Eq. © , the decay width T is determined self-consistently 
by r/2 = — Im [II(a>)]. The self-consistent calculation be- 
gins by putting L/2 = 70 into the right-hand side of 
Eq. ([8]). By summing the right-hand side, we have a new 
r/2 via r/2 = — Im [n(w)] and we put it into the right- 
hand side again. This calculation is repeated until II(w) 
is converged. We use Eq. © with J7 = g a su(uj)/£ for Vl 
in Eq. ([8]) . Here g a ff is the off-site electron-phonon matrix 
element and u(ui) is the amplitude of the LO mode. We 
adopt g ff = 6.4 eV4^ A similar value is obtained by a 
first-principles calculation with the local density approxi- 
mation. 46 We use a harmonic oscillator model which gives 
u(ui) — yfh/2M c N u LO where M c is the mass of a carbon 
atom and N u is the number of hexagonal unit cells. 

In Fig. [3]Ja), we plot the renormalized energy Kloq + 
Re [II(u;o)] as a function of Ep for the LO and TO modes 
in a N = 9 Z-NR at room temperature (300K). Since 
the TO mode decouples from electron-hole pairs, the 
self-energy IT(ct;o) vanishes and the frequency of the TO 
mode does not change from ujq. On the other hand, the 
LO mode exhibits a KA effect. The error-bars extend- 
ing up to ±r/2 in Fig. [3] represent the broadening of 
the phonon frequency due to the finite life time of the 
phonon. The LO mode shows the largest broadening 
effect of r « 5cm -1 when Ep — OeV. The value of L 
decreases quickly when Ep > O.leV. This is because of 
the Pauli exclusion principle by which a resonant decay 
of the LO mode is forbidden when Ep > 

For comparison, we show the renormalized energies of 
the LO and TO modes for a (5, 5) armchair SWNT as 
a function of Ep by the black curves in Fig. ©a). The 
TO mode exhibits no broadening but a softening with 
a constant energy ~ —30 cm -1 . The absence of broad- 
ening is due to the fact that the Bloch function can be 
taken as a real number for lowest energy sub-bands, i.e., 



for vertical transition denoted by the dashed line 1 in 
Fig. ©b) even when a Z-NR is rolled to form an arm- 
chair SWNT. The details are given in Sec. |lVl For the 
LO mode, the broadening of the Z-NR is smaller than 
that of the armchair SWNT because \Vp\ of a Z-NR is 
smaller than that of armchair SWNTs for the lower en- 
ergy bands (see Fig. (He)). In fact, the real part of the 
right-hand side of Eq. © is a negative (positive) value 
when E e h > Uuiq (E e h < Thoq). Thus, electron-hole pairs 
with higher (lower) energy contribute to the frequency 
softening (hardening) , 40 i 47 Since the energies of the edge 
states for Z-NR are smaller than Tiloq, the edge states 
may contribute a frequency hardening like the one shown 
around Ep = for a (5, 5) SWNT. The absence of the 
hardening confirms that the edge states do not contribute 
to n(w ) because is negligible. 
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FIG. 3: (color online) (a) The _Bf dependence of the renor- 
malized energies of the LO and TO modes in a N — 9 Z-NR 
(red) and of the LO and TO modes in a (5, 5) SWNT (black), 
(b) The -Bf dependence of the renormalized frequencies of the 
LO and TO modes in a N = 9 A-NR (red) and of the LO and 
TO modes in a (9,0) SWNT (black). 

It is noted that, for the renormalized energies of the 
LO and TO modes in NRs shown in Fig. [31 we have not 
included all the possible intermediate electron-hole pair 
states created by a given phonon mode in evaluating ^ ch 
in Eq. ©. For example, vertical transition denoted by 
the arrow 3 in Fig. Ufa) may be included as a possible 
intermediate state in evaluating X) h m ^q- © although 
such intermediate state does not satisfy the momentum 
conservation. In this paper, we do not consider the con- 
tribution of momentum non-conserving electron-hole pair 
creation processes in evaluating ^ oh in Eq. ©. 



B. Armchair NRs 

Next we study the KA in A- NRs. The zigzag SWNTs 
are cut along their axis and flattened out to make A-NRs. 
It is known that one third of zigzag SWNTs exhibit a 
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metallic band structured It is interesting to note that 
if we cut the bonds along the axis of a metallic zigzag 
SWNT in order to make an A-NR, the obtained A-NR 
has an energy gap. Namely, unrolling a metallic (3z,0) 
SWNT results in a N = 3i — 1 A-NR with an energy gap. 
Instead, unrolling a semiconducting (3i + 1,0) SWNT 
results in a gap-less N = 3i A-NR and unrolling a semi- 
conducting (3i+2, 0) SWNT results in a N = 3i + l A-NR 
with an energy gap. The one-third periodicity of metal- 
licity is maintained even if zigzag SWNTs are unrolled by 
cutting the bonds. Since metallicity is a necessary con- 
dition for the KA, we study the KA in N = 3i metallic 
A-NRs here. 

In order to specify the lattice structure of an A-NR, 
we use integers I € [0, N] and J £ [0, M] in Fig.ffla). In 
a box specified by IJ in Fig. 0Ja), there are two A atoms 
and two B atoms. For convenience, we divide the two 
A (B) atoms into up-A (up-B) and down-A (down-B), 
as shown in Fig. [Ha). The wave function then has four 
components: (e'( l2< ' J /\/M)(^ A , <^ B , <^ A , <^ B )' where 
k is the wave vector along the armchair edge. 
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FIG. 4: (a) The lattice structure of an A-NR. (b,c) The 
eigenvectors of the Y point LO and TO phonon modes are 
illustrated. The LO mode satisfies <57o,i = —257o,2 and 
<57o,2 = 570,3. The TO mode is characterized by 570,1 = 
and 570,2 = —570,3- 

In Fig. [3Jb) and (c), we show phonon eigenvectors 
of the r point LO and TO modes, respectively. The 
eigenvector of the LO (TO) mode is parallel (perpen- 
dicular) to the armchair edge. The LO mode satisfies 
^70,1 = —2^7o, 2 and #70, 2 — <57o,3, while the TO mode 
satisfies ^70,1 = and £70, 2 = — ^70,3- Since £70,2 and 
^70,3 are perturbations that do not mix <p u and 4>d, the 
electron-hole pair creation matrix element can be divided 
into the following two parts: 



7 A V ( + <570,2<5/',/+l 

h.c. 



1,1' 



y dB 



h.C. 



(9) 



where the Hermite conjugate (h.c.) of V u (V d ) is defined 



as (570,3(5/',/ + <57o, 2<5/',/- 
can rewrite Eq. ((9| as 



((57o,2(5/',/ + (57o, 3 <5/',/+i). We 



Va = 2i^2 (*70, 3 Im [^IKb] + <ho,2lm [^VSb]) , 
/ 

V d = 2*]T (^7o, 2 Im [0 d 10 dB ] + <57o, 3 Im ) . 



(10) 



The perturbation <57o,i mixes (f)^ and 4>^ as 
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/ 0uA \ f 
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ik2e 5^ ,x 
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(11) 



so that U can be rewritten as 



U = i2(5 7o , 



E 

/ 



(im [< 



,ik2iil* 



Im 



1 1* 



(12) 



Now, it can be shown that each energy eigenstate sat- 
isfies the following equations (see Appendix [A]) . 



/ 



/ 



/ 



^dA^dB; 



E^ 



1+1 

dA^dB • 



(13) 



Due to these conditions, the TO mode causes a special 
cancellation between V u and Vd as V n + Vd — since 
<57o,2 + (57o,3 = 0. In addition, we obtain U = from 
<57o,i = 0. Thus the TO modes in A-NRs decouple from 
electron-hole pairs and do not undergo a KA. For the LO 
mode, on the other hand, there is no cancellation between 
Vu and V d , and the matrix element for the LO mode is 
given by V L = U + V u + V d . By setting £70, 2 = £70,3, 
—2(57o,3 = £7 and <57o,i = (57, we calculate Vl and put it 
into Eq. to obtain U(ujq). 

The solid curves in Fig. [3£b) give the E-p dependence 
of the renormalized frequencies hu>o + Re [II(a>o)] for the 
LO and TO modes in a N — 9 A-NR at room temper- 
ature. The frequency of the TO mode is given by ojq 
showing that the TO mode decouples from electron-hole 
pairs. The LO mode undergoes a KA. For comparison, 
we show the renormalized frequency of the LO and TO 
modes in a (9, 0) zigzag SWNT as the black curves in 
Fig. [3Jb). It is found that |n(o;o)| for the LO mode of 
a N = 9 A-NR is smaller than |JI(cc?o)| for the LO mode 
of a (9, 0) zigzag SWNT. It is because there are two lin- 
ear energy bands near the K and K' points in metal- 
lic zigzag SWNTs, while there is only one linear energy 
band in metallic A-NRs, and the KAs in A-NRs are sup- 
pressed slightly as compared to those in metallic zigzag 
SWNTs. We note that the broadening in A-NRs is still 
larger than that in Z-NRs because of the absence of the 
edge states near the armchair edges. The TO mode of 
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a (9,0) zigzag SWNT is down shifted. However, there 
is no dependence on E-p, which indicates that only high 
energy electron-hole pairs contribute to the self-energy. 
We will explore the KA effect for the TO mode in zigzag 
SWNTs in Sec. US 

We have considered NRs with a long length (10/xm) 
in calculating the self-energies H(wq) shown in Fig. [21a) 
and (b). For NRs with short lengths, the effect of the 
level spacing on H(u>) is not negligible. For example, the 
level spacing in armchair SWNTs becomes 0.12eV when 
L = 30nm, which is comparable to Tlujq/2. Thus the level 
spacing affects the resonant decay. The effect of the level 
spacing on the KAs in NRs will be reported elsewhere. 

III. RAMAN INTENSITY 
A. Raman Activity 

In the Raman process, an incident photon excites an 
electron in the valence energy band into a state in the 
conduction energy band. Then the photo-excited elec- 
tron emits or absorbs a phonon. The matrix element 
for the emission or absorption of a phonon is given by 
the el-ph matrix element for the scattering between an 
electron state in the conduction energy band and a state 
in the conduction energy band, which is in contrast to 
that for the el-ph matrix element for electron-hole pair 
creation which is relevant to the matrix element from a 
state in the valence energy band to a state in the conduc- 
tion energy band. This matrix element for the emission 
or absorption of a phonon is given by removing the mi- 
nus sign from — </>jg (or — </>^ B dB ) of the final state in the 
electron- hole pair creation matrix elements in Eqs. fl}, 
©, and (fTTj) . This operation is equivalent to replacing 
Im (Re) with Rc (Im) in Eqs. ([5]). Q, (fTU |) and flT g) . 

As a result, the Raman intensity of the TO (LO) modes 
in Z-NRs is proportional to Re [0 a *0b] (Im [0a* ^b] )■ 
Thus the TO modes are Raman active, while the LO 
modes are not. Because the TO modes in Z-NRs are free 
from the KA, the G band Raman spectra exhibit the 
original frequency of the TO modes, Hluq. For A-NRs, 
on the other hand, the cancellation between V u and Va 
occurs for the TO modes, and the TO modes are then 
not Raman active, while the LO modes are Raman ac- 
tive. Since the LO modes in A-NRs undergo KAs, the 
renormalized frequencies, Kloq + II(a>o), w iU appear be- 
low the original frequencies of the LO modes by about 
20 cm -1 (see Fig. [3]). Thus, the G band spectra in A- 
NRs can appear below those in Z-NRs, which is useful 
in identifying the orientation of the edge of NRs by G 
band Raman spectroscopy. Our results are summarized 
in TABLE [J combined with the results for armchair and 
zigzag SWNTs. 

For the Raman intensity of armchair SWNTs, we ob- 
tain the same conclusion as that of Z-NRs, that is, the 
TO modes are Raman active, while the LO modes are 
not. For metallic zigzag SWNTs, on the other hand, the 



cancellation between V u and Vd which occurs for the TO 
modes in A-NRs is not valid. Then the TO modes, in 
addition to the LO modes, can be Raman active. How- 
ever, as we will show in Sec. lIVl since the matrix element 
for the emission or absorption of the TO modes vanishes 
at the van Hove singularities of the electronic sub-bands, 
then we can conclude that the TO modes are hardly ex- 
cited in resonant Raman spectroscopy. It is interesting to 
compare these results with another theoretical results for 
the Raman intensities of SWNTs. In Ref. 48, it is shown 
using bond polarization theory that the Raman inten- 
sity is chirality dependent. In particular, for an armchair 
(zigzag) SWNT, the A lg TO (LO) mode is a Raman ac- 
tive mode, while the A\ g LO (TO) mode is not Raman 
active. These results for SWNTs are consistent with our 
results. We think that it is natural that the Raman in- 
tensity does not change by unrolling the SWNT since the 
Raman intensity is proportional to the number of carbon 
atoms in the unit cell and is not sensitive to the small 
fraction of carbon atoms at the boundary. 

TABLE I: Dependences of the KAs and Raman intensities 
on the r point optical phonon modes in NRs and metallic 
SWNTs. O an< A x represent 'occurrence' and 'absence', re- 
spectively. A means that the KA is possible, but the broad- 
ening effect weakens due to the presence of the edge states. 
V means that the KA is possible, but the -Ef dependence is 
suppressed by the decoupling from the metallic energy band 
crossing at the Dirac point. 



mode Kohn anomaly Raman active 



Z-NRs 


LO 


A 


X 




TO 


X 


O 


A-NRs 


LO 


o 


o 




TO 


X 


X 


Armchair SWNTs 


LO 


o 


X 


(rolled Z-NRs) 


TO 


V 


o 


Zigzag SWNTs 


LO 


O 


o 


(rolled A-NRs) 


TO 


V 


X 



B. Comparison With Experiment 

We prepare graphene samples by means of the cleav- 
age method to observe the frequency of the G band 
Raman spectra. In many cases, graphene samples ob- 
tained by the cleavage method show that the angles be- 
tween the edges have an average value equal to multiples 
of 30°. This is consistent with the results by You et 
al— Figure El[a) shows an optical image of the exfoliated 
graphene with the edges crossing each other with an an- 
gle of ~ 30°. This angle can be considered as evidence of 
the presence of edges composed predominantly of zigzag 
or armchair edges. Note that, the obtained sample is 
characterized as a multi-layer graphene. We estimated 
the number of layers to be about five based on the be- 
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havior of the G' band. The sample is placed on a SiC>2 
(100) surface 300 nm in thickness. 

The Raman study was performed using a Jobin-Yvon 
T64000 Raman system. The laser energy is 2.41eV 
(514. 5nm), the laser power is below lmW and the laser 
spot is about 1 /im in diameter. Figure [5](b) shows spec- 
tra for finding the position dependence of the G band. 
The results show that the G band frequency depends 
on the position of laser spot. When the spot is focused 
near the upper edge (A) or far from the edge, the posi- 
tion of the G band is almost similar to that of graphite 
(1582cm" 1 ). However, a softening of the G band is 
clearly seen when the laser spot moves to the vicinity 
of the lower edge (B). 

Based on our theoretical studies, we find that the G 
band observed near the upper edge consists only of the 
TO mode, while the G band observed near the lower edge 
comes from the LO mode, because the G band near the 
lower edge shows a down shift which is considered to be 
due to the KA effect of the LO mode. Then, we speculate 
that the upper (lower) edge is dominated by the zigzag 
(armchair) edge. 

It should be noted that our experiment does not prove 
that the observed down shift of the G band near the 
lower edge is due to the KA effect, since we do not ex- 
amine the E-p dependence. There is a possibility that 
the observed downshift of the G band is related to me- 
chanical effects. Mohiuddin et alM- observed that the 
G band splits into two peaks due to uniaxial strain and 
both peaks exhibit redshifts with increasing strain. The 
edge in this work somehow has half of it suspended and 
this may decrease the vibration energy. This effect may 
explain why for the upper edge in Fig. [5] the frequen- 
cies are slightly downshifted compared with the spectra 
taken at the center of the graphene sample. Moreover, it 
is probable that the physical edge in this work is a mix- 
ture of armchair and zigzag edges ^ Note also that we 
measured the D band in order to confirm that the iden- 
tification of the orientation of the edge is consistent with 
the fact that the D band intensity is stronger near the 
armchair edge than the zigzag edge i 25 ' 27 i 41 However, we 
could not resolve the difference of the intensity near the 
upper and lower edges clearly. Zhou et al^S observed by 
means of high-resolution angle-resolved photo-emission 
spectroscopy experiment for epitaxial graphene that the 
D band gives rise to a kink structure in the electron self- 
energy and pointed out that an interplay between the 
el-ph and electron-electron interactions plays an impor- 
tant role in the physics relating to the D band. We will 
consider these issues further in the future. 



IV. PSEUDOSPIN 

In the preceding sections we have shown for Z-NRs that 
the LO modes are not Raman active and that the TO 
modes do not undergo KAs. These results originate from 
the fact that the Bloch function is a real number. Besides, 



(a) 




(b) 




1550 1580 1610 
Raman shift [cm -1 ] 



FIG. 5: (color online) (a) An optical image of a graphene 
sample. The sample is characterized as a multi-layer graphene 
(~ 5 layers) sample, (b) The position dependence of the G 
band frequency. From the theoretical results obtained in this 
paper, we conclude that the G band appearing near the upper 
(lower) side of the edge consists only of the TO (LO) mode, 
and that the upper (lower) edge is dominated by the zigzag 
(armchair) edge. 



the TO modes in A-NRs are not Raman active and do 
not undergo KAs. This is due to the cancellation between 
V u and Vd for the TO modes. The LO modes in A- 
NRs undergo KAs since the Bloch function is a complex 
number. The absence or presence of a relative phase 
between the <f>\ and </>b of the Bloch function is essential 
in deriving our results. In this section, we explain the 
phase of the Bloch function in terms of the pseudospin, 51 
and clarify the effect of the zigzag and armchair edges on 
the phase of the Bloch function. 
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A. Absence of a Pseudospin Phase in Z-NRs 

Here we use the effective-mass model in order to un- 
derstand the reason why the Bloch function in Z-NRs is 
a real number. In the effective-mass model, the Bloch 
functions in the conduction energy band near the K and 
K' points are given by-Si 



<PK(k x , k y ) = — 



(p K '(k x , k y ) 




(14) 



where k x and k y (k' x and k') are measured from the K 
(K') point and the angle 6 (8') is defined by k x + ik y = 
\k\e i9 (k' x + ik y = |k'|e ie '). It is noted that k x (k y ) is 
taken as parallel to the zigzag (armchair) edge. Then k y 
is reflected into —k y at the zigzag edge, and the scattered 
state is given by 



y>K(kx, -k y ) 



(15) 



as shown in Fig. [6l 




zigzag edge 



armchair edge 



FIG. 6: In fc-space, we consider a state near the K point (solid 
circle) and states which are scattered by the zigzag and arm- 
chair edges (empty circles). k x (k y ) is taken as parallel to 
the zigzag (armchair) edge. The arrows in the insets indi- 
cate the direction of the pseudospins. The (a y ) component of 
the pseudospin is reversed at the zigzag edge, while the (a x ) 
component is preserved at the armchair edge. 

The relative phase of the wave function at the A and 
B sublattices can be characterized by the direction of 
the pseudospin. The pseudospin is defined by the ex- 
pectation value of the Pauli matrices <J x ,y,z with re- 
spect to the Bloch function.— For <fiK{k x ,k y ), we have 
(<r x ) = cos9, (a y ) = si) 
we have (a x ) — cos 8, ( 



{a z ) = 0. For ip K (k x , - H y ) , 
= -sin0, (<r z ) = 0. Then 



due to an interference between the incoming and reflected 
waves, we have (<r y ) = for the Bloch function near the 
zigzag edge. The condition (a y ) = means that the 
Bloch function becomes a real number. In fact, since 
K> = 2£jlm[?^] and (a x ) = 2£ J Re[ ( ^] 
for the Bloch function ((/r£, (j>^) in Eq. (O, the el-ph ma- 
trix elements Vt and Vl (Eq. ©) are proportional to 
(a y ) and (<r x ), respectively. In Appendix IBl we show 
the relationship between the Bloch function c^k in the 
effective-mass model (Eq. (1141) ) and the Bloch function 
(</>a,</>b) in the tight-binding model (Eq. 0). 

We point out that the condition (a y ) = is not satis- 
fied in the case of armchair SWNTs (rolled Z-NRs) ex- 
cept for the lowest energy sub-bands of k y = 0. This is 
the reason why we see in Fig. El^a) that the TO mode 
exhibits no broadening but a softening with a constant 
energy ~ —30 cm -1 in a (5,5) SWNT. It is interesting 
to note that the Aharanov-Bohm flux applied along the 
tube axis shifts the cutting line a 52 ' 53 i 54 and can make 
k y for the lowest energy sub-bands nonzero. Thus the 
Aharanov-Bohm flux makes that (a y ) — sin 8 ^ even 
for the lowest energy sub-bands, for which the TO mode 
can exhibit a broadening^ 

Since the effective-mass model describes the physics 
well in the long wave length limit, an advantage in the 
above discussion of using the effective-mass model is that 
it is not necessary for the edge of graphene NRs to be well 
defined on an atomic scale in order that we have a can- 
cellation of (cr y ). This may be a reason why we observe a 
softening of the G band near the edge of an armchair-rich 
sample experimentally as shown in Sec. IIII Bl 



B. Coherence of the Pseudospin in A-NRs 

On the other hand, a state near the K point with 
(k x , k y ) is reflected by the armchair edge into a state near 
the K' point with (k' x , k' y ) — (—k x ,k y ), and the scattered 
state is given by (pK'(—k x , k y ). In this case, by putting 
8' = 7r — 8 into tfK'(k X) k y ) in Eq. (jT5J) , we obtain 



<PK'{-kx, k y ) 




(16) 



(<j y ) flips at the zigzag edge as shown in Fig. Thus, 



which is identical to the initial Bloch function, 
iPK.(k x ,k y ). Thus the relative phase between the A and 
B Bloch functions is conserved thorough the reflection 
by the armchair edge so that the Bloch function can not 
be reduced to a real number. Namely, the reflections by 
the armchair edge preserve the pseudospin as shown in 
Fig. [6] It is expected that the relative phase makes it 
possible that KAs occur not only for the LO mode but 
also for the TO mode near the armchair edge. However, 
as we have shown in Eq. (|13[) . the armchair edge gives 
rise to the cancellation between V u and Vd, so that the 
TO modes in A-NRs do not undergo KAs. We consider 
whether Eq. (fT3"|) is satisfied in the case of zigzag SWNTs 
or not, in order to see if the KA effect is present in zigzag 



9 



SWNTs or not. By applying the Bloch theorem to zigzag 
SWNTs, we have 



<PLb = (e 4(/+1/2)fe 72V^B, 

= ( e l(/+1/2)fc 72ViV)^ A , 



(17) 



where we set ip^, = * (<£>a, Vb)- Using these equations, we 
obtain 



J 

'AdA'PdB = ^ e 2 VA^B. 

J 



(18) 



Thus the first equation in Eq. (fT?|) is not satisfied for 
zigzag SWNTs. Similarly, we have 



1 



-e 2 Va^b, 



E^dA^B 1 = J^Va^E, 



(19) 



which shows that the second equation in Eq. (|13[) is not 
fulfilled, either. Thus the TO modes can undergo KAs 
because the cancellation between V u and Vd is not possi- 
ble for zigzag SWNTs. In fact, by putting Eqs. (fT8|) and 
(TT9"j) into Eq. (TTUD, we S et 



K + Vd = «(<57o,2 + £70,3) sin 9 cos 
+ *(^7o,3 - ^70,2) cos 9 sin 



ka 

T 

ka 



(20) 



where we set ip\(fiB = e 10 ■ Because the TO modes satisfy 
^70,2 + ^70,3 = 0, V u + Vd can take a nonzero value for 



iu ,- Kl = /2f - )-,|.;; Sill ( ^ I C"sf 



(21) 



It is noted that Eq. (|2~lj) vanishes when 9 = ±ir/2. This 
condition 9 — ±tt/2 is satisfied for low energy electron- 
hole pairs when the energy band crosses the Dirac point. 
In other words, high energy electron-hole pairs in the 
sub-bands can contribute to a frequency softening of the 
TO mode in zigzag SWNTs. This is why we obtain the 
down shift of the TO mode for a (9,0) zigzag SWNT 
as shown in Fig. [3J In "metallic" zigzag SWNTs, the 
curvature effect shifts the position of the cutting line 56 
of the metallic energy band from the Dirac point and 
produces a small energy gap . 57 i 58 ' 59 i 60 In this case, the 
low energy electron-hole pairs satisfy cos 9 ^ in Eq. (|2"Tj) 
and they contribute to a frequency hardening of the TO 
modes in metallic zigzag SWNTs. The curvature effect 
gives rise to a change of the Fermi surface and results in 
KAs for the TO modest 



Similarly, the matrix element for the emission or ab- 
sorption of the TO modes in zigzag SWNTs is given by 



sin ( 



(22) 



which does not vanish in general. This shows that the TO 
modes in zigzag SWNTs can be Raman active. However, 
since the van Hove singularities of sub-bands in zigzag 
SWNTs are located on the k x axis (and satisfy = 0), 
the sin# factor tells us that the TO modes are hardly 
excited in resonant Raman spectroscopy. 



V. DISCUSSION 

The theoretical analysis performed in this paper is 
based on the use of a simple tight-binding method which 
includes only the first nearest-neighbor hopping inte- 
gral and its variation due to the atomic displacements. 
The approximation used is partly justified because the 
deformation-potential and the el-ph matrix element with 
respect to the second nearest-neighbors is about one or- 
der of magnitude smaller than that of the first nearest 
neighborsj 45 ' 46 However, we have not considered the ef- 
fect of the overlap integral. The overlap integral breaks 
the particle-hole symmetry and may invalidate our re- 
sults. Besides, we have neglected the contribution of mo- 
mentum non-conserving electron-hole pair creation pro- 
cesses in evaluating Yleh m Eq- ©• Although this is an 
approximation which works well for thin NRs, the inclu- 
sion of the momentum non-conserving electron-hole pair 
creation processes may invalidate our results. We will 
elaborate on this idea in the future. 



VI. SUMMARY 

In summary, the LO modes undergo KAs in graphenc 
NRs while the TO modes do not. This conclusion does 
not depend on the orientation of the edge. In Z-NRs, 
the Raman intensities of the LO modes are strongly sup- 
pressed because the wave function is a real number, and 
only the TO modes are Raman active. As a result, the 
KA for the LO mode in Z-NRs would be difficult to ob- 
serve in Raman spectroscopy. In A-NRs, only the LO 
modes are Raman active owing to the cancellation be- 
tween V u and Vd- The "chirality" dependent Raman in- 
tensity derived for NRs is the same as the chirality depen- 
dent Raman intensity for SWNTs calculated in Ref. HI. 
The strong down shift of the LO mode makes it possi- 
ble to identify the orientation of edges of graphene by 
the G band Raman spectroscopy due to the "chirality" 
dependent Raman intensity. 
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APPENDIX A: DERIVATION OF EQ. (fTS)) 



By rewriting the Bloch function of Z-NRs ^ A , 
Eq. as 



2i \g 



C{g,<t>) 
2i 



e(s>0)- 



(Bl) 



a iJ<j> 
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where c.c. denotes the complex conjugation of the first 
term, one sees that (4>\, 0b) is a real number as a result 
of the cancellation of the imaginary part between the first 
and second terms. By introducing a new Bloch function 

(p<f, as 



In this Appendix, we derive Eq. (j 1 3|) by using mir- 
ror and time-reversal symmetries. Let us introduce the 
mirror-reflection operator M by 



(° 












«a\ 








1 







€b 





1 










Ka 


\1 















(I = 0,...,N), (Al) 



where k is the wave vector along the armchair edge and 
ii is the band index. By applying M to the energy 
eigen-equation H k 4> n , k = E ntk (f> n> k, we get MH k (f> n , k = 
E n ,kM(j) n ^ k . Since the Hamiltonian satisfies MH k M~ x = 
H^ k , we obtain M(f> n ^ k = e l ^cf> ny - k where is a phase 
factor. 

On the other hand, due to the time-reversal symmetry, 
we have 0* fc = e 1 ^ (j> n ,-k- Thus, by combing the time- 
reversal symmetry (0* k = e*^ (j> n .-k) with the mirror 
symmetry (M</>„, fc = '4> n -k) , we get 



M<j) nk 



(A2) 



that is, 



(Kb\ 




(Kl\ 








Ka 


= e^" 


Kb 


Kb 




Kl 


Wua) 




\KbJ 



(A3) 



Using this condition, one sees that Eq. (fl"3|) is satisfied 



APPENDIX B: RELATIONSHIP BETWEEN 
EQ. © AND EQ. (fUD 

Here we will show for Z-NRs that the Bloch function 
derived using the tight-binding lattice model (Eq. ([7])) is 
a superposition of incoming and reflected Bloch functions 
derived using the effective-mass model (Eq. (JHJl). 



= °{9,4>) c *fg + e 

2i 9 I £ (5» 



the Bloch function (0 A , 0g) is expressed by 



(B2) 



<p<pe 



-iJ<fi 



(B3) 



Because of the different signs in the exponents of e %J ^ 
and e~ tJ ^ in Eq. (|B3|) . 4> may be thought of as the wave 
vector perpendicular to the zigzag edge (k y ) multiplied 
by a lattice constant (£) as <f> — k y £. Assuming that 
<j> = k y £, the zigzag edge reflects a state with k y (<^e iJ< ^) 
into a state with — k y (ip- < pe~ lJ ^). We then expect that 
the Bloch function near the zigzag edge is given by 



p^ + p.^ 1 ^. 



(B4) 



It is noted that Eq. (|B4j) is different from Eq. (jB3|) be- 
cause f—(f, is not identical to p*^ in general. However, we 
will get ip-ct, = ip^ for Z-NRs because we may assume 
that the normalization constant C(g, (f) in Eq. (|B2p sat- 
isfies C(g,—4>) = —C*(g,(/>) without loss of generality. 
Therefore, Eq. (1B4|) is consistent to Eq. (|B3j) . which in- 
dicates that the assumption (cj> — k v £) is appropriate. 
The condition p-^, — <p\ is a non-trivial condition since 
it is satisfied only for the zigzag edge and is essential for 
(</>A, 0b) ^° ^ e a rea ^ number. 

Using Eq. (fl4|) we obtain ip-&{k x , —k y ) — p>K(k x ,k y )* 
which corresponds to ip—cj, = ip\. In fact, by putting 
(.9 + e~ # ) = \g + e^|e~ ie into Eq. £52} and by using 



t(9,4>) 
as 



\9- 



44 



the Bloch state p><p can be written 



V<t> 



V2 (e« 



(B5) 



where O and 9 are the same as each other near the K 
point. 



* Email address: SASAKI. Kenichi@nim s.go.jp1 



1 X. Li, X. Wang, L. Zhang, S. Lee, and H. Dai, Science 319, 



11 



1229 (2008). 

2 L. Jiao, L. Zhang, X. Wang, G. Diankov, and H. Dai, Na- 
ture 458, 877 (2009). 

3 D. V. Kosynkin, A. L. Higginbotham, A. Sinitskii, J. R. 
Lomeda, A. Dimiev, B. K. Price, and J. M. Tour, Nature 
458, 872 (2009). 

4 R. Saito, M. Fujita, G. Dresselhaus, and M. S. Dresselhaus, 
Appl. Phys. Lett. 60, 2204 (1992). 

5 Z. Liu, K. Suenaga, P. J. F. Harris, and S. Iijima, Phys. 
Rev. Lett. 102, 015501 (2009). 

6 X. Jia, M. Hofmann, V. Meunier, B. G. Sumpter, 
J. Campos-Delgado, J. M. Romo-Herrera, H. Son, Y.-P. 
Hsieh, A. Reina, J. Kong, et al, Science 323, 1701 (2009). 

7 C. O. Girit, J. C. Meyer, R. Erni, M. D. Rossell, 
C. Kisielowski, L. Yang, C.-H. Park, M. F. Crommie, M. L. 
Cohen, S. G. Louie, et al., Science 323, 1705 (2009). 

8 K. Kobayashi, Phys. Rev. B 48, 1757 (1993). 

9 M. Fujita, K. Wakabayashi, K. Nakada, and K. Kusakabe, 
J. Phys. Soc. Jpn. 65, 1920 (1996). 

10 K. Nakada, M. Fujita, G. Dresselhaus, and M. S. Dressel- 
haus, Phys. Rev. B 54, 17954 (1996). 

11 K. Wakabayashi, M. Fujita, H. Ajiki, and M. Sigrist, Phys. 
Rev. B 59, 8271 (1999). 

12 J. M. Holden, P. Zhou, X.-X. Bi, P. C. Eklund, S. Bandow, 
R. A. Jishi, K. D. Chowdhury, G. Dresselhaus, and M. S. 
Dresselhaus, Chem. Phys. Lett. 220, 186 (1994). 

13 A. M. Rao, E. Richter, S. Bandow, B. Chase, P. C. Eklund, 
K. W. Williams, S. Fang, K. R. Subbaswamy, M. Menon, 
A. Thess, et al., Science 275, 187 (1997). 

14 M. Sugano, A. Kasuya, K. Tohji, Y. Saito, and Y. Nishina, 
Chem. Phys. Lett. 292, 575 (1998). 

15 M. A. Pimenta, A. Marucci, S. A. Empedocles, M. G. 
Bawendi, E. B. Hanlon, A. M. Rao, P. C. Eklund, R. E. 
Smalley, G. Dresselhaus, and M. S. Dresselhaus, Phys. 
Rev. B 58, R16016 (1998). 

16 Z. Yu and L. Brus, J. Phys. Chem. B 105, 1123 (2001). 

17 A. Jorio, R. Saito, J. H. Hafner, C. M. Lieber, M. Hunter, 
T. McClure, G. Dresselhaus, and M. S. Dresselhaus, Phys. 
Rev. Lett. 86, 1118 (2001). 

18 S. Bachilo, M. Strano, C. Kittrell, R. Hauge, R. Smalley, 
and R. Weisman, Science 298, 2361 (2002). 

19 M. Strano, S. Doom, E. Haroz, C. Kittrell, R. Hauge, and 
R. Smalley, Nano Letters 3, 1091 (2003). 

20 S. Doom, D. Heller, P. Barone, M. Usrey, and M. Strano, 
Appl. Phys. A: Mater. Sci. Process. 78, 1147 (2004). 

21 M. S. Dresselhaus, G. Dresselhaus, R. Saito, and A. Jorio, 
Physics Reports 409, 47 (2005). 

22 L. G. Cangado, M. A. Pimenta, B. R. A. Neves, 
G. Medeiros-Ribeiro, T. Enoki, Y. Kobayashi, K. Takai, 
K.-i. Fukui, M. S. Dresselhaus, R. Saito, et al., Phys. Rev. 
Lett. 93, 47403 (2004). 

23 A. C. Ferrari, J. C. Meyer, V. Scardaci, C. Casiraghi, 
M. Lazzeri, F. Mauri, S. Piscanec, D. Jiang, K. S. 
Novoselov, S. Roth, et al., Phys. Rev. Lett. 97, 187401 
(2006). 

24 J. Yan, Y. Zhang, P. Kim, and A. Pinczuk, Phy. Rev. Lett. 
98, 166802 (2007). 

25 M. A. Pimenta, G. Dresselhaus, M. S. Dresselhaus, L. G. 
Cancado, A. Jorio, and R. Saito, Phys. Chem. Chem. Phys. 
9, 1276 (2007). 

26 C. Casiraghi, A. Hartschuh, H. Qian, S. Piscanec, 
C. Georgi, A. Fasoli, K. S. Novoselov, D. M. Basko, and 
A. C. Ferrari, Nano Letters 9, 1433 (2009). 

27 L. Malard, M. Pimenta, G. Dresselhaus, and M. Dressel- 



haus, Physics Reports 473, 51 (2009), ISSN 0370-1573. 

28 H. Farhat, H. Son, G. G. Samsonidze, S. Reich, M. S. 
Dresselhaus, and J. Kong, Phys. Rev. Lett. 99, 145506 
(2007). 

29 K. T. Nguyen, A. Gaur, and M. Shim, Phys. Rev. Lett. 

98, 145504 (2007). 

30 Y. Wu, J. Maultzsch, E. Knoesel, B. Chandra, M. Huang, 
M. Y. Sfeir, L. E. Brus, J. Hone, and T. F. Heinz, Phys. 
Rev. Lett. 99, 027402 (2007). 

31 A. Das, A. K. Sood, A. Govindaraj, A. M. Saitta, 
M. Lazzeri, F. Mauri, and C. N. R. Rao, Phys. Rev. Lett. 

99, 136803 (2007). 

32 W. Kohn, Phys. Rev. Lett. 2, 393 (1959). 

33 O. Dubay, G. Kresse, and H. Kuzmany, Phys. Rev. Lett. 
88, 235506 (2002). 

34 S. Piscanec, M. Lazzeri, F. Mauri, A. C. Ferrari, and 
J. Robertson, Phys. Rev. Lett. 93, 185503 (2004). 

35 M. Lazzeri and F. Mauri, Phys. Rev. Lett. 97, 266407 

(2006) . 

36 S. Piscanec, M. Lazzeri, J. Robertson, A. C. Ferrari, and 

F. Mauri, Phys. Rev. B 75, 035427 (2007). 

37 N. Caudal, A. M. Saitta, M. Lazzeri, and F. Mauri, Phy. 
Rev. B 75, 115423 (2007). 

38 S. Pisana, M. Lazzeri, C. Casiraghi, K. S. Novoselov, A. K. 
Geim, A. C. Ferrari, and F. Mauri, Nature Materials 6, 198 

(2007) . 

39 T. Ando, J. Phys. Soc. Jpn. 77, 014707 (2008). 

40 K. Sasaki, R. Saito, G. Dresselhaus, M. S. Dresselhaus, 
H. Farhat, and J. Kong, Phys. Rev. B 77, 245441 (2008). 

41 Y. You, Z. Ni, T. Yu, and Z. Shen, Applied Physics Letters 
93, 163112 (2008). 

42 R. Saito, G. Dresselhaus, and M. Dresselhaus, Physical 
Properties of Carbon Nanotubes (Imperial College Press, 
London, 1998). 

43 K. Sasaki, S. Murakami, and R. Saito, J. Phys. Soc. Jpn. 
75, 074713 (2006). 

44 K. Sasaki, S. Murakami, R. Saito, and Y. Kawazoe, Phys. 
Rev. B 71, 195401 (2005). 

45 J. Jiang, R. Saito, G. G. Samsonidze, S. G. Chou, A. Jorio, 

G. Dresselhaus, and M. S. Dresselhaus, Phys. Rev. B 72, 
235408 (2005). 

46 D. Porezag, T. Frauenheim, T. Kohler, G. Seifert, and 
R. Kaschner, Phys. Rev. B 51, 12947 (1995). 

47 K. Sasaki, R. Saito, G. Dresselhaus, M. S. Dresselhaus, 

H. Farhat, and J. Kong, Phys. Rev. B 78, 235405 (2008). 

48 R. Saito, A. Jorio, J. H. Hafner, C. M. Lieber, M. Hunter, 
T. McClure, G. Dresselhaus, and M. S. Dresselhaus, Phys. 
Rev. B 64, 085312 (2001). 

49 T. M. G. Mohiuddin, A. Lombardo, R. R. Nair, A. Bonetti, 
G. Savini, R. Jalil, N. Bonini, D. M. Basko, C. Galiotis, 
N. Marzari, et al., Phys. Rev. B 79, 205433 (2009). 

50 S. Y. Zhou, D. A. Siegel, A. V. Fedorov, and A. Lanzara, 
Phys. Rev. B 78, 193404 (2008). 

51 K. Sasaki and R. Saito, Prog. Theor. Phys. Suppl. 176, 
253 (2008). 

52 H. Ajiki and T. Ando, J. Phys. Soc. Jpn. 62, 2470 (1993). 

53 S. Zaric, G. N. Ostojic, J. Kono, J. Shaver, V. C. Moore, 
M. S. Strano, R. H. Hauge, R. E. Smalley, and X. Wei, 
Science 304, 1129 (2004). 

54 E. D. Minot, Y. Yaish, V. Sazonova, and P. L. McEuen, 
Nature 428, 536 (2004). 

55 K. Ishikawa and T. Ando, J. Phys. Soc. Jpn. 75, 084713 
(2006). 

56 R. Saito, K. Sato, Y. Oyama, J. Jiang, G. G. Samsonidze, 



12 



G. Dresselhaus, and M. S. Dresselhaus, Phys. Rev. B 72, 
153413 (2005). 

R. Saito, M. Fujita, G. Dresselhaus, and M. S. Dresselhaus, 
Phys. Rev. B 46, 1804 (1992). 

C. L. Kane and E. J. Mele, Phys. Rev. Lett. 78, 1932 



(1997). 

T. Ando, J. Phys. Soc. Jpn. 69, 1757 (2000). 

L. Yang and J. Han, Phys. Rev. Lett. 85, 154 (2000). 



